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We reconsider post-inflation gravitino production, in the context of hidden sector super- 
gravity models. We discuss the possible role of supersymmetry breaking from finite tem- 
perature effects in enhancing the rate for gravitino production and argue that there is 
no such enhancement. Our conclusion is based on a simple decoupling argument, which 
is independent of temperature. We also characterize the properties of goldstino-like hy- 
drodynamic fluctuations that arise in supersymmetric models at finite temperature. We 
show that they cannot lead to enhanced gravitino emission via infrared divergences. We 
comment on an analogy with axion emission from hadronic matter. 
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1. Introduction 

Supersymmetry is an attractive theoretical framework as it provides a solution to 
the hierarchy problem for fundamental scalar masses. With dynamical mechanisms for 
supersymmetry breaking, the gauge hierarchy problem (mw "C mpi) can also be explained. 

It is of great interest then to study the effects of supersymmetry on various experi- 
mental and cosmological data. There is of course as yet no direct evidence for the existence 
of the supersymmetric partners of the particles of the Standard Model. We must then rely 
on studying indirect effects. There has been much work done in the last two decades in 
these regards. In this paper, we will focus on certain cosmological issues. 

In particular, the standard cosmological models predict that the early Universe was 
filled with a plasma; the constituents of this plasma have densities determined by equilib- 
rium thermodynamics. Since the Bose-Einstein and Fermi-Dirac distributions are different, 
the plasma is populated by different amounts of on-shell fermions and bosons. 

In a theory in which supersymmetry plays a role,@ temperature effects may invalidate 
various cancellations implied by the symmetry between fermions and bosons. It is impor- 
tant then to investigate what role supersymmetry really plays in the early universe, and 
to study possible effects of temperature on various calculations in supersymmetric models. 

In all extant models of supersymmetry breaking, the theory consists of two sectors. 
In the first sector, supersymmetry is broken by some mechanism which in general we will 
not specify. The observable particles, i.e., those of the standard model, are in a second 
sector. The two sectors are separate in the sense that they interact only very weakly with 
one another. In the standard hidden sector scenarios, this interaction is provided by non- 
renormalizable gravitational effects. In models of low-energy supersymmetry breaking, 
such as those recently discussed by Dine and Nelson[|l], the communication of the two 
sectors is mediated by gauge interactions. We note that in the latter models, the gravitino 
can be extremely light (m 3 / 2 = F/mpi > ra^/mpi), and cosmology places no bounds on 
its interactions, as long as it is sufficiently light. We will thus concentrate in this paper 
on hidden sector models. 

We have been particularly motivated by the recent work of W. Fischler[]3[ on post- 
inflation gravitino production, in which it is suggested that temperature effects can greatly 
enhance the Goldstino component of gravitino production cross-sections, leading to very 
tight constraints. The plan of this paper is as follows. In the next section, we consider 

1 Here we mean that supersymmetry is unbroken at some scale in the zero-temperature theory. 
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the equivalence theorem for Goldstinoes and gravitinoes at high energies/temperatures 
and give a decoupling argument based on finite-temperature Ward identities that suggests 
that Goldstino amplitudes will in fact not be enhanced. In the third section, we suggest 
a calculation that can correctly compute the total Goldstino production rate, making 
proper use of real-time finite-temperature quantities. This calculation avoids any direct 
discussion of individual Feynman diagrams; we argue that the total rate can be computed 
in (re-summed) perturbation theory for weak coupling. In the second half of this paper, we 
will review older literature on the effects of temperature on super symmetric theories. We 
clarify and extend results on how supersymmetry is broken, and its physical manifestation 
as propagating modes in the plasma. 

2. Goldstino-Gravitino Equivalence and Decoupling 

In this section we consider the possible relevance of gravitino production in the pri- 
mordial plasma after reheating. Since we are dealing with temperatures T > m 3 / 2 , the 
transverse and longitudinal modes of the gravitino field may play very different roles. In 
particular the transverse (spin 3/2) component couples always with gauge strength 1/mpi. 
Thus the rate of change of its density n 3 / 2 is at most of order T 6 /m Pl . In conventional 
cosmological scenarios, this small rate gives the usual bounds. If the gravitino is stable, it 
must be lighter than a few keV, in order that it not dominate the energy density. ||J] If it 
does decay, in the absence of dilution by inflation, nucleosynthesis will constrain m 3 / 2 ^ 10 
TeV.[|J In inflationary scenarios, gravitinoes will be diluted, with 713/2/% ~ TpH/ m Pi- 
However if the reheat temperature Trh after inflation is too high (Trh > 10 9-10 GeV), 
gravitinoes will re-attain sufficient density as to have problematic effects on nucleosynthe- 
sis, entropy release after nucleosynthesis, light element photo-dissociation and distortions 
in the microwave background. || 

The longitudinal (spin 1/2) mode, however, can interact more strongly. Indeed the 
interactions of this mode at E > m 3 / 2 , are well described by those of the Goldstino 
via the equivalence theorem. At high energies they are proportional to a power of 
(E/mpim 3 /2) ~ E/F, and could perhaps || lead to a large n 3 / 2 oc T 8 /F 2 . The cos- 
mological consequences could then be rather strong. In what follows, however, we argue 
that in the elementary particle models of phenomenological interest the rate of goldstino 
production is never so enhanced. Indeed it is always suppressed by at least one power of 
1/M 2 , where M is the scale of the physics communicating SUSY breaking to conventional 
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matter. In particular, when SUSY breaking is communicated by gravity, we have M ~ rripi 
and the rate for production of helicity 1/2 gravitinos is never more important than the 
rate for helicity 3/2. 

The models of interest consist of a hidden and an observable sector. It is assumed 
that in the absence of interactions between them the hidden sector breaks supersymmetry 
while the observable one does not. The fundamental Lagrangian of the models of interest 
may thus always be written in the form: 



£ = CoM) + C hid (X) + £e n £ n (0,X) (2.1) 



where e is a small parameter. In models where the interaction between the two sectors 
is induced by some physics at a scale M we have e = 1/M. For instance in conventional 
supergravity models M = mp\. The dynamics of Chid breaks SUSY at a scale F = M%. 
The SUSY breaking effects in the observable sector are then 0(F/M), and it must be 
F/M ~ Mz for phenomenological reasons. It is manifest from the form of this Lagrangian 
that, in the limit e — > 0, the Goldstone fermion of spontaneously broken supersymmetry 



does not couple to the fields </>. The form of eq. (2^1) then shows explicitly that, as long 
as perturbation theory is valid, also at finite temperature the rate of Goldstino production 
from observable matter is suppressed by at least e 2 . We note that the Goldstino production 
cross-section of Ref. || does not satisfy this property. We also note that if we were interested 
in the production of Goldstinos from a thermal bath of hidden sector particles we would 
indeed get a rate scaling like T 8 /F 2 . However this is not a situation of cosmo logical 
interest. 

By performing a local field redefinition, the Goldstino coupling with matter can be 
written in terms of the supercurrent. At lowest order in the Goldstino x field we have 



(2.2) 



where in an obvious notation S£ are the contributions to the supercurrent from the terms 
in eq. (^.1|). The decoupling of goldstino and observable sector is not explicit in eq. 
(|2.2|) . When calculating scattering amplitudes the decoupling shows up via cancellations 
among diagrams which are related to each other by supersymmetry. The diagrammatic 
way of reasoning is however very dangerous when trying to draw conclusions on finite 
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temperature processes. Indeed the presence of a thermal bath breaksa supersymmetry 
(as well as Lorentz invariance), and one might then naively expect that the "effective" 
diagrams will be related to each other up to terms proportional to a power of T. That 
is, since the static limit of propagators are affected by finite temperature, it is tempting 
to say that the mass splittings should be replaced by temperature dependent factors. As 
a result of this, the suppression of processes with goldstinos and ordinary matter would 
be 0(T/F) rather then e ~ in clear contradiction with our decoupling argument. 

However, when using eq. (|2.2|) to parametrize the goldstino couplings, one should note 
that is a conserved current in the limit e — > 0. This means that the decoupling, in 
this language, is the result of a Ward identity. Real time correlation functions will obey 
Ward identities that are identical to their zero temperature counterparts, except that the 
correlation functions should be interpreted as finite temperature averages. In particular, 
the current is indeed conserved in the limit M — > oo, and so any correlation function 
involving a Goldstino and observable matter must also vanish in this limit, irrespective 
of temperature. The introduction of ensemble parameters like temperature (or chemical 
potential for abelian symmetries), corresponds to the choice of a particular state, which 
may not be invariant under some symmetries of the system. Ward identities, on the other 
hand, express properties of the dynamics, i.e., the Hamiltonian, and as such do not depend 
on the particular state or ensemble average. In other words, the introduction of ensemble 
parameters can lead to a Goldstone representation of a symmetry but not to its explicit 
breakdown. Then an effective finite temperature description of real time processes, if any 
is possible, must preserve the form of the Ward identities. This is reminiscent of what 
happens with chiral symmetry in QCD, where the low energy pion lagrangian realizes the 
same Ward identities of the fundamental QCD dynamics. 

Now one should note that such decoupling arguments can in general be violated by 
the effects of resonant enhancement. This occurs when a relatively large mixing arises 
amongst nearly degenerate states, because of a small energy denominator. However, at 
finite temperature, because of damping phenomena, coherence can be maintained only for 
finite times < tdamp, and the resonant enhancement is limited. Indeed the only modes 
that are undamped are those corresponding to Goldstone singularities at zero energy and 
momentum. It is known that finite temperature Ward identities (even in the absence of 
spontaneous supersymmetry breaking) require the existence of a fermionic excitation with 

2 See section 4. 



4 



some properties similar to a Goldstino. In the next two sections we will show that these 
Goldstino-like hydrodynamic modes cannot lead to enhanced gravitino production, simply 
because of the form of their dispersion relation. In the Appendix, we comment on the 
case of axion emissioni from thermal pions, where a partial enhancement does indeed take 
place due to resonant mixing with a Goldstone mode. 

3. Goldstino production at finite temperature: a calculation 

We now present the outline of a calculation of the production rate of gravitinoes at 
finite temperature. It is well known that the interpretation of finite temperature Green's 
functions as corresponding to scattering amplitudes is tricky. Essentially the problem is 
one of boundary conditions: we do not have simply asymptotic states to deal with, but a 
full plasma. The plasma is assumed to consist of thermalized visible sector particles; the 
hidden sector must be (nearly) unpopulated. Otherwise, the Goldstino couples directly to 
those particles and will thermalize with that sector, leading to the usual gravitino problem. 
This is clearly a condition on inflation models: the inflaton, whose decay is responsible for 
reheating the universe, must not couple to the hidden sector (except through unavoidable 
gravitational couplings). This is sufficient to ensure that the decay of the inflaton will 
thermalize the visible sector fields, but not the hidden sector. 

In such a situation, we wish to study the approach of the Goldstino to thermal equi- 
librium. The idea that we will use here is to attempt to compute the inclusive production 
rate,@ where we sum over all final states including a Goldstinoi and average over initial 
states by using a thermal density matrix pp = e~@ H /Z. We thereby avoid having to 
interpret individual scattering amplitudes. 

The Goldstino couples to the visible sector through the supercurrent, with 



where the coefficients a, b are of order one. The total transition probability per unit 4- 
volume for production of single gravitinoes may be written: 






(3.2) 



1 



/ 



d 3 q 



I 



d 4 x e iq - x tr 4 (O(x)O(0)) p . 



M 2 



(2tt) 3 2£ 9 



We thank S. Thomas and W. Fischler for bringing this example to our attention. 

4 We thank T. Banks for discussions on this point. 

5 States with more than one Goldstino will be parametrically suppressed. 
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The use of the r.h.s. of eq. ( |3.1|) is justified here, since we are bracketing on on-shell 
states, which satisfy the full equations of motion. In the first line of ( |3.2| ) , we have written 
the result in terms of Dirac spinors and is the Goldstino 4-momentum. Since we are 
interested in T ^> m 3/ / 2 we have taken q 2 = 0. We note that the thermal matrix element 
(O(x)O(0)) p = Tr [p^O(x)(5(0)] is a function of real times, but is not a time-ordered 
quantity. Thus we can also write this expression in terms of the discontinuity (or the 
imaginary part) of the gravitino self energy 

w x = / —4^ 5^ tr 6 21m (O0) q a = (3-3) 

x M 2 J (2n) 3 2E q e?" + 1 % V /q ' P M 2 K ! 

We are thus lead to an expression for the total production rate which is suppressed 
by 1/M 2 as long as the matrix element is non-singular in the limit M — > oo, i.e., as 
long as there is no resonant mixing. In what follows we will argue that this is indeed the 
case. Notice that a behaviour w x ~ (1/M) r with < r < 2 would be already physically 
interesting; the result of Ref. corresponds to r = 0. 

In the limit M — > oo, the (thermalized) visible sector, in which (OO)p is calculated, 
corresponds to an exactly supersymmetric theory: the MSSM without soft supersymme- 
try breaking terms. The singularity required to modify the 1/M 2 behaviour of eq. ( |3.3| ), 
would be physically interpreted as a mode in a fermionic channel O, which propagates 
undamped on the light-cone. The only modes which propagate undamped at finite T, 
according to prime principles, are Goldstones corresponding to spontaneously broken sym- 
metries. In fact these have to be undamped at (u,k) = (0,0), in order to saturate the 
broken Ward identities. For (a;, k) ^ (0, 0) we have no physical reason to expect an un- 
damped mode in an interacting theory, and in fact experience shows that these modes 
are always damped. It is known that finite temperature leads to a spontaneously broken 
representation of supersymmetry and the appearance of a Goldstino-like mode. Indeed 
fermionic Green's functions like (OO)p are singular at (u,k) = (0,0). Fortunately it is 
possible to characterize in a quite general way the properties of this singularity. As shown 
in the next section, its dispersion relation in a weakly coupled theory at high T is given 
by uj/k = 1/3, i.e., strictly off the light-cone. In other words, in the limit in which the 
Goldstino is undamped, only the real part of (OO) (and not the imaginary one) has a 
1/k singularity when approaching (u,k) = (0,0) from the light-cone. (Indeed, cfr. (3.3), 
a much stronger singularity would be needed to affect w x ). We will thus conclude that 
no resonant mixing can arise with the zero-temperature Goldstino in eq. ( |3.3j ), so that 
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w x ~ 1/M 2 . In the Appendix we discuss an example where the Goldstone singularity gets 
arbitrarily close to the light-cone as k — > 0, so that a deviation from the naive behaviour 
for a total emission rate does indeed take place. 

Indeed, since w x is well defined as M — > oo it will also be calculable in perturbation 
theory at weak (gauge and Yukawa) coupling. The leading contributions come from two- 
loop diagrams involving A-terms (cf. the first term in eq. ( |3.1| )), as well as contributions 
from the gaugino mass (second term in eq. (|3.1|) ). One-loop diagrams do not contribute to 
w x in the susy limit: they correspond to two body decay processes with goldstino emission, 
which are forbidden by phase space. The differential rate d 3 w x /d 3 /c, when computed at 
momenta < gT requires some care to be taken. The matrix element (|3.3|) will in general 
contain infrared divergences when calculated in perturbation theory. These divergences 
signal a breakdown in naive perturbation theory when the energy carried by a given line 
(or by all lines entering a vertex) in a graph is of order g s T; fortunately, it is known how to 
resum the perturbation series. H The result, for nonzero external momenta (u,q) is that 
all modes are screened at order g s T, including gauge particles. This is consistent with 
our general argument that w x is finite. We stress once more that in order to avoid the 
1/M 2 suppression for Goldstino production, w x (even after resummation) would have 
to be infinite (we have taken M — > oo). There is one remaining subtlety in the case of 
non-Abelian gauge theory; namely in the integration region u, q < g 2 T, we expect further 
infrared divergences. These are similar to those encountered in the static limit, where for 
example, the calculation of the free energy of hot QCD irrevocably breaks down at some 
order.! Thus ultimately, we cannot calculate the partial rate in this regime. If something 
dramatic were to happen to the rate in this momentum region, it would happen in a 
perturbatively incalculable way. But again, consistent with the general physical argument 
discussed above, we feel that this is implausible. For instance, it is known that in hot 
QCDfllOl] the quasiparticle spectrum in the quark propagator can lead to enhancements 
in the partial rate at low energies. However, it seems reasonable that higher order effects 
lead to a finite width for these modes which would cut off any zero-energy singularity and, 
if we may use the analogy here, w x would then indeed be finite. 

6 For recent progress on this problem, see Ref. Q. 
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4. Temperature Goldstinoes 

In this final section, we will discuss some aspects of the behaviour of supersymmetric 
theories at finite temperature and the appearance of thermal Goldstino modes. Supersym- 
metry is broken at finite temperature because of the non-zero energy density in the ground 
state; the breaking thus appears spontaneous and there exists a fermionic Goldstone mode. 
In this section we study the physics of this mode; we begin by considering theories in which 
supersymmetry is unbroken at zero temperature. Much of what we have to say is a review 
of portions of the literature, although we believe some of our comments are new. 



In Ref. it is argued that temperature treats fermions and bosons differently, so 
supersymmetry is expected to be broken. They suggest that it is explicitly broken by the 
boundary conditions (the antiperiodic boundary condition for fermions leads to a problem 
in imaginary time Ward identities: the Grassmann variational parameter e cannot be taken 
to be a constant) and that no massless fermion mode need arise. 

In Refs. |T^],[T^], it was realized that in order to study dynamical effects, one must 
consider the problem in the real-time finite temperature formalism. It was established 
that the Ward identities at finite temperature are essentially those of the zero temperature 
theory with matrix elements interpreted as finite temperature correlation functions. In 
particular, Boyanovsky ||12|| established a Ward identity involving a supercurrent- fermion 
matrix element and found a fermionic pole with no mass gap when (F)p is non-zero. 

Let us consider the problem now in more general terms. The superalgebra relation 



{<Q a ,Ql} 



implies that a state is invariant under supersymmetry if and only if its energy is zero. 
Then, when the ground state itself has non-zero energy, supersymmetry is spontaneously 
broken and a Goldstone fermion is present in the spectrum. The existence of the Goldstino 
is established by writing eq. ( |4.1| ) in its local form, as a Ward identity 

d,(T (S£(X)S£{0))) = * 4 (*) (OTK* (4-2) 

where <S M and Q^ v are respectively the supercurrent and the stress energy tensor. By 
integrating eq. ( |4.2| ) , a non-zero righthand side implies the existence of a zero momentum 
singularity in the current-current correlator. Moreover, Lorentz invariance of the vacuum 
constrains the singularity to be of the form 1/ft, signaling the presence of a massless 
fermion. Notice, though, that eq. (|4.2|), is a local operator identity determined by the 



dynamics of the system, and, as such, is valid on any stationary, translation invariant 
state. In particular in a thermal bath, (9) ^ 0, and we expect long range correlations in 
the system, and the absence of a mass gap. The absence of Lorentz invariance, though, 
does not allow one to establish the nature of the singularity in a model independent way 
ll~4| . For instance in the free massless Wess-Zumino model, the left hand side in eq. (|4.2| ) 



is easily evaluated and the singular part has the same form as the Klimov-Weldon [|15[ 
self-energy from hard thermal loops. The form is involved, with logarithmic singularities 
on the light-cone, but the only "pole-type" singularity is at the single point u = 0, k = 
(where u and k are respectively energy and three- momentum) . Thus there is long-range 
order, but without a propagating wave. 



In fact, as also noticed in Ref. [|16|j , eq. ( [I.2|) is the supersymmetric analogue of the 
Ward identity expressing the spontaneous breakdown of Lorentz invariance in a thermal 
bath 

d A x dx(T(M Xtl „(x)6 pv (y))}{3 = (g PP uri (y) - g vp Q m {y) + g m Q pv {y) - g^O P p(y)) p (4.3) 

where M^p is the Lorentz current. In a thermal bath, (9 PV ) qk. g prj and the right-hand side 
is non-zero. Then the physical interpretation of Eq. (|4.3|) is that there are propagating 
hydrodynamic fluctuations (sound waves) corresponding to variations in the stress-energy 
density. These modes correspond to local boosts of the thermal bath, which in the limit of 
infinite wave-length cost zero energy.S In addition we also know from hydrodynamics that 
these modes correspond to propagating waves only at wave-lengths that are much bigger 
than the mean free path (see for example ref. |17|]). In particular, in a free field theory, 
the mean free path is infinite, and eq. ( fL3| ) does not give rise to any interesting physical 
wave at finite (u, k). This is consistent with our findings for eq. (|4.2|) . By analogy, it also 
suggests that in an interacting supersymmetric theory a hydrodynamic fermionic wave will 
describe the fluctuations of the system above some long, but finite, wavelength. It also 
suggests, that, by a reasoning that parallels the one in hydrodynamics, one should be able 
to derive the wave equation for the Temperature Goldstino in a quite general way. 

Let us briefly recall the procedure in hydrodynamics. One considers a fluid that 
in stationary conditions in its restframe is described by a stress energy tensor Q pv = 
pg^v + (p + p)go/j,goi>, where p and p are respectively the pressure and energy density. 



7 Notice the subtlety in counting Goldstone modes for this space-time symmetry; only one of 
three modes propagates. 
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The Goldstone mode of broken Lorentz (or Galilean) symmetry corresponds to small local 
boosts of the quiescent system which are described by a space-time dependent velocity 
three- vector Vi. To lowest order in v, the energy momentum in the "global" rest frame of 
the bath will have the boosted form 

Onu = <V + A v [<V] = P9tiu + (p + p) (gov9ov - Vig itI go„ - Vig 0fX g iu ) (4.4) 

The continuity equation for 9 together with the equation of state of the fluid, then give 
the equation of propagation of the Goldstone mode,i the sound velocity equaling v s = 
y^dp/dp^ 

Let us now consider the case of supersymmetry. In a thermal bath the density of 
supercurrent (S£) vanishes at equilibrium. In the presence of a space time dependent 
Goldstino oscillation £ the supercurrent S is given by 

3£ = {Qa, §£} = e^(a^(x)) & (4.5) 

and conservation of the current corresponds to the Goldstino wave equation. We thus 
obtain the dispersion relation 

u = (-) k. (4.6) 



.9 , 

The velocity of the Goldstino vq = p/p should be compared to the velocity of sound 
v s = \Jdp/dp. In a weakly coupled theory we have respectively in the non-relativistic and 
relativistic regimes: vq = = T/m and Vq = Vg = 1/3. In the intermediate regimes, 
though, the relation vq = does not hold. 

Eq. (|4.6|) agrees with the explicit calculation of Ref. [0, where the fermion self- 



energy in the Wess-Zumino model was studied at low temperature. However we suspect 



8 Notice that to derive the equation the derivatives of p, p have to be considered non zero and of 
order v, i.e., the longitudinal part of the Goldstone field v "mixes" with other scalar modes. This 
is because when Lorentz invariance is broken, there can be mixing terms between the Goldstone 
field and heavier modes which contain just one and not two time derivatives. For example, in the 
liquid helium case, the Lagrangian contains the term p^dtn where p is the chemical potential, 
$ the real part of the condensate and tt the Goldstone phase. Integrating out the massive $ 
produces a relevant quadratic piece for n. As we will see this phenomenon does not take place for 
Goldstinoes due to the linearity in (dt,d x ) of the wave equation. 

9 Here we are considering zero chemical potential, so that sound waves are associated with 
fluctuations in the energy density. 
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that in the treatment in Ref. [18], and also [12], [13], higher order corrections are not 
completely under control. Though these effects will not change eq. (|4.6|) (apart from 
the form of higher order corrections to the vacuum energy and pressure) , we suspect that 
they can significantly affect the residue of the Goldstino pole. In what follows we briefly 
motivate our doubts. For the sake of the argument we will stick to the low temperature 
case which was studied in Ref. jl8],|12],|13|]. Eq. ( [4.2|) suggests that the order parameter 
for supersymmetry breaking is a composite operator, i.e., the stress-energy tensor. We 
thus expect the Goldstino to lie predominantly in a composite channel, suggesting that 
the description of this mode in terms of the elementary fields in the lagrangian is indeed 
non-perturbative. In Refs. Jl2|] , Jl3[] , it was noticed that in the massive Wess-Zumino 
model, (0) is displaced from its zero temperature value (po = m/g and we have a vacuum 
expectation value for the auxiliary field (F) = d^W ^ generated at finite temperature; 
thus the Goldstino has a non-zero overlap with the elementary fermion field in the theory. 
The fermion and boson masses, as calculated from the tree lagrangian, are thus split 
by A = m A — = — m 2 B = gF, where A, B and if) are respectively the scalar, 
pseudoscalar and fermion. The 1-loop effective fermion self energy T,(uj,k) is consistently 
calculated by using the free propagators with the masses corrected by the background in- 
field value. At zero momentum in real time one finds S ~ 5 (2 /A. In our case, we must have 
m + E(0, 0) = as required by the broken-supersymmetry Ward identity. [0 The dominant 
contributions to E(a>, k) for (u, k) <C (A/T, A/vmT), are determined by denominators of 
the form 

- - (4 7) 

El - (E A -u) 2 A + 2pk - 2E A u V ' ; 

where p is the integration 3-momentum. The above denominator is the same as that of 
hard thermal loops, apart from the "infrared regulator" A. It is however clear that any 
higher order contribution affecting the fermion and boson propagators can affect the result, 
since it is determined in lowest order by the "small" A. For instance, the genuine 1-loop 
corrections to the scalar masses at low temperature (i.e., not those coming from (F) ^ 0) 
are not smaller than order A. Of course the fact that individual corrections affect the result 
does not mean that the full correction will, especially since SUSY cancellations might be 
expected. Nonetheless it is not clear to us how to implement the "nominally" higher 
order effects in a manifestly SUSY consistent manner. This is especially clear at high 
temperature where the 1-loop corrections to propagators dominate the tree level. When 
one is interested in calculating processes at low momenta uj ~ gT the resummation of hard 
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thermal loops [19] is the correct approach. However, for ultralow momenta 10 <C g 2 T we 
are not aware of any developed technique. 

The damping rate of this mode was also calculated in Ref. |T3[ and was found to be 



exponentially small at small momentum^ 7^ ~ exp(—/3A/2k), where A is the effective 
fermion-boson mass splitting induced by temperature. This exponential suppression arises 
presumably only at one loop. At this order the attenuation of the Goldstino is determined 
by a three body process, in which a boson (fermion) in the bath absorbs a thermal Goldstino 
and becomes a fermion (boson). Since fermions and boson effective masses have a splitting 
A, energy-momentum conservation requires an energy 0(A/uj) for the absorbing particle, 
thus the Boltzmann suppression at low 10. At higher orders the thermal Goldstino can be 
absorbed in a four- (or higher) body process. In this case, we do not expect the above 
phase space exponential suppression, and presumably the only suppression at low (u, k) 
will be due to the fact that the Goldstino is derivatively coupled. 

The dispersion relation was obtained in Ref. |l8j by expanding around (u, k) = 0, and 
is apparently not sensitive to A. Fortunately, for the discussion of gravitino production, 
we need only have control of the dispersion relation ( |4.6| ). The fact that the "singularity" 
lies far off the light cone is clear evidence that it cannot play any relevant role in gravitino 
production, as the mode will not mix with the gravitino. Given this, the behaviour of the 
Goldstino damping 7^ cannot play any role in our discussion. Nonetheless it would be 
physically interesting to find a way to calculate this quantity. 

Our analysis should be contrasted with the case studied in the Appendix, where 
the mode that leads to enhanced axion production is given by the pion. It that case, 
indeed, the relevant mode is not a purely collective one but a thermally "dressed" massless 
particle. Thus, not surprisingly, the mass shell is very close to the light-cone and a partial 
enhancement takes place. 

We would like to thank T. Banks, M. Dine, G. Farrar, R. Pisarski, S. Raby, and S. 
Thomas for useful discussions. We are indebted to the Aspen Center for Physics where 
this work was begun. 



10 Of course we know that the damping rate must be exactly zero at cj, k = 0, since there must 
be a Goldstone singularity. 
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5. Appendix 

In this Appendix we show the details of the pion-axion mixing example which was 
mentioned in the text. We focus on the case of two light quark flavors u, d. 

The coupling between the axion field a and hadronic matter is given, at lowest order, 
by the (a, tiq) mass matrix 

^ = to) (™r ™r) ) = B(a, tt ) ( m /i m , m ^ d ) ( I ) (5.1) 

where f a ^> are the axion and pion decay constants, while B =< qq >~ Aq Cd . We 
are interested in the rate for axion production from a thermal bath of hadrons. To lowest 
order in m 2 n and using the formalism of Section 3, the rate of axion production per unit 
time and volume is 

dn a 4 f d 3 k a I 




At - m " 



I (2^ki^H7rFT !Disc[G " (£ "' fc " )] (5 ' 2) 



Where G^n is the pion propagator in the thermal bath. This function has been calculated 
at finite temperature up to 0(1/ f^) in the chiral expansion. |f20f Even in the chiral limit, 



both the real and imaginary parts of G~^.(uj, k) are non- vanishing on the light-cone. The 
only surviving singularity is at (u, k) = (0,0), consistent with Goldstone's theorem. The 
thermal self-energy H n7T (uj, k) can be written as 

U^iu;, k) = (R + il) ^k 2 + ^0(K 2 : ml) + i^O(K 2 , ml) (5.3) 

J TT J TT J TT 

where K 2 = uj 2 — k 2 , while R and / are numerical coefficients and R > 0. In order to 
illustrate the infrared enhancement effect we consider the limit <C T <C f^. In this 
regime the dominant contribution comes from relativistic pions. Then, as will become 
clear below, we need only consider the term proportional to R + il in eq. ( |5.3|) . The other 
terms are subleading with respect to the tree-level chirality breaking terms in G^. The 
discontinuity of the pion propagator then reads 

Of 4 Tb 2 

iDiscG^(E a , k a ) = 3-2 - 2 (5.4) 



The above expression behaves like 1/k 2 in the range k a 3> fci n f ~ m. K (f^/T yR). Then, 
for k a in this range, the integral in eq. ( |5.2| ) is linearly infrared divergent, related to the 
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bosonic nature of the pion, and its proximity to the light-cone. This divergence is cut-off 
at momenta of order k- m f. Neglecting numerical factors, the total rate behaves like 



dn o f 2 _ 

— ocm„ -ocms (5.5) 

at m^l 

displaying an enhanced and non-analytic behaviour in the quark mass m. Notice also the 
growth of eq. ( |5.5|) at low temperature. As temperature is lowered the pion density gets 
smaller, but the overlap of a pion excitation with a massless particle improves due to the 
smaller real and imaginary parts in IT,,-,,-. Then there is coherence on a longer scale and a 
pion is more likely to oscillate into an axion. 

Of course, the above results can be qualitatively interpreted by considering the quan- 
tum mechanics of pion-axion oscillations taking into account damping phenomena. For 
a pion emitted at t = 0, the probability of finding an axion at time t is given by 
P(t) — e 2 s'm 2 (Am 2 t/2E), where e is the mixing angle and eAm 2 = m^. However, in 
the presence of damping, this will be limited by the finite coherence length of the pions. 
For a damping rate F 3> m 2 /E, we have P^^a ~ m„ T / '(EF) 2 . Then we simply find the 
rate of axion production 

dUa mL Tn^ (5.6) 



dtd 3 k E 2 T 2 

qualitatively consistent with the above more formal analysis. 
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